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Abstract. Gravity is nearly a universal constant in the cusp of an NFW galaxy halo. Inside 
this external field an isothermal gas sphere will collapse and trigger a starburst if above a 
critical central pressure. Thus formed spheroidal stellar systems have Sersic-profile and satisfy 
the Faber- Jackson relation. The process is consistent with observed starbursts. We also recover 
the Mbh — a* relation, if the gas collapse is regulated or resisted by the feedback from radiation 
from the central BH. 
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1. Tight correlation of formation of Black Hole and Bulges 

The formation of central black holes (BHs) in galaxies is likely a rapid process since most 
quasars have already formed at redshift z > 2. The co-relation between the BH mass and the 
velocity dispersion of the spheroidal (bulge) component is so tight that it is hard to explain 
unless bulges form as fast as BHs to keep their growth neck-to-neck. While at the present day 
the BH accretion rate is completely decoupled from the bulge growth, it is possible that their 
growth was correlated during the violent feedbacks at high redshift. Indeed starburst activities 
peak at similar redshifts as the quasars as a whole. 

In two recent papers by Xu, Wu, & Zhao (2007) and Xu & Wu (2007), we propose that bulges 
can form by a rapid collapse due to radial instability of isothermal gas. This model has the nice 
feature of forming bulges before disks. 

Here we iterate the key steps of the above scenario, but without invoking the gravothermal 
instability as in Xu et al. (2007). Instead we follow the example of Elemgreen (1999) and find 
the equilibrium configurations of the maximum gas mass inside an external gravity. We also 
generalize the argument to a star-gas mixture to show that we can form bulges with a reasonable 
profile. Assuming the rapid star bursts in the bulge are regulated by the accretion-driven wind 
of the central black hole, we derive the BH mass-stellar dispersion relation. 



2. A universal constant gravity scale for dark halos 

In the Cold Dark Matter (CDM) framework, baryons fall into the potential well of CDM, cool 
and condense into stars. Here we consider the properties of gas and stellar equilibrium in the 
external field of dark matter. 

The background dark matter distribution is often described by the NFW density distribution 
for dark matter (Navarro, Frenk & White 1997), which has a density pnfw ~ p s r s /r inside a 
scale radius r s . In the central region which concerns the galaxy bulge, we note an interesting 
universal scale for the dark halo gravity 

jdm(i") = — — = 27rGTI ~ 10 msec a, a ~ 1 (2.1) 
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where n = p s r s ~ 130Mqpc _2 H is a column density and S(M V i r , z, c) ~ 1 is a shallow function of 
the halo virial mass M v i r , the redshift, and the concentration c. In another words, the gas where 
formed the bulge was imbedded in an uniform external field from the dark matter potential. We 
can also define a dark matter central pressure 

P DM = g 2 DM /(AnG) = goMll/2 = p DM (r) ■ gnu ■ (r/2) (2.2) 

for later use. 



3. Maximum gas mass sustainable by halo gravity 

In general for a gas and a stellar sphere imbedded in an external DM gravity, Consider 
imbedded in an external DM gravity an isothermal gas sphere p(r) of sound speed a, and an 
isotropic stellar sphere p* (r) of dispersion <j* in quasi-static equilibrium, the potential at a given 
radius r is 

, f r ( GM + GM, \ 

Jo \ 9DM + H ) (3 } 

The equilibrium satisfies the equations 

9»m + - 2 ] = Tr = -a. , (3.2) 



dM(r) _ dM,{r) 



p(r)dr p*(r)dr' 



(3.3) 



where we define a\ = (l + r)<r 2 . Here a position-independent feedback factor T 2> 1 is introduced 
because the radiation from the star burst and the accreting central black hole can generate an 
additional opacity-induced pressure (pcr 2 )r on the dusty gas sphere (but not the stellar sphere), 
countering the gravity. 

First consider the stage where the stellar mass is negligible, so M* <C M. Rewrite the equations 
in term of the following dimensionless mass, radius and density, 

x(m) = — , p(m)= 2 H \/ 1 (3.4) 



crfgDu/G' a\ /gnu' g 2 DM (4irG)- 

and express the rescaled gas mass m as the independent coordinate, the problem is recasted to 
solving the pair of dimensionless ODEs 

x 2 dp(m) m x 2 dx(m) 1 

dm x(m) 2 ' dm p(ni) 

For each value of p(0), the gas density profile under the hydrostatic equilibrium can be totally 
determined with the following initial conditions at the center for the radius x(0) =0 and the 
rescaled density 

2 2 

„m\ P° ai d — 9dm /o a \ 

po = p(0) = — , Pdm = — ^, (3.6) 

Fdm 47rG 

where p(0)/(l + F) equals the ratio of the gas central pressure po"" 2 vs the dark matter's pressure 
under self-gravity Pdm- 

Computing the gas equilibrium for a range of core pressure po (see figure [TJ, we find that 
the gas density generally falls montonically with radius or mass. All models have finite mass 
out to infinite radius where p = 0. This interesting behavior is due to the deep potential well 
of the external gravity, which makes the isothermal density drop exponentially with radius as 

~Tdm ^ n P p M ^ ~ r (-^0i hence the mass converges quickly if neglecting self-gravity. 

The finite mass of these gas spheres will give another interesting behavior. There is a critical 
core pressure 

Po = f^«30 (3.7) 



How to form bulges/ ellipticals in dark halos as fast as central black holes? 



3 



above which the gas density p(M) of a parcell of gas dM no longer increases monotonically with 
an increase of central pressure, and in fact the total mass will decrease with increasing po after 
it reaches a maximum value 

M max « 4.3 ( — ^-^) . (3.8) 

These limits on gas central pressure and total mass are related to the instability first discussed 
by Elmgreen (1999). Gas sphere above certain critical mass Mm ax oc a\ or critical central gas 
density or pressure do not have stable solutions; adding tiny amount of gas would lead to collapse. 
It is interesting to speculate that the bulge formation originates from such a gas instability. 



4. Post-starburst Mass Profile 

Our models are also generalizable while gas is converting to stars. The simplest solution of 
eqs. (3.2-3.3) would be a model where stars trace gas radial distribution, so we have 

P_M = M 1 (r 1 = _ k _ aj = 

p(r) M{r) 1 - /. ' a 2 + ^ ' 1 ' 

where the position-independent factor /*(£) is the fraction of gas formed into stars at time t. 
Such a solution is possible if the feedback is regulated by star formation. Rescaling the gas-only 
solution (eq. 3.7), we obtain the stability criterion 

P0<j2 " P *' 0(7 * < 3QP DM , P DM = (4.2) 



M? = ~ 4 x 10 11 M s 9 _ - , (4.3) 

qdmG \200kms 1 ' 



(i-/.)/(i + r) f, ^ " m AtvG- 

Gas can turn into stars quasi-staticly where maintaining the above equality at the critical 
density and mass. Such formed stellar system would have a total mass 

„ = i^_„ 4 xlP 11 " / 
9dmG 

and a profile 

M (r) r y r 
*\ L = / j(y)(Wdy), y=^-^, j(y) « OMy- 1 exp(-1.6 2/ 1/1 - 2 ), (4.4) 

where j(y) is numerically fitted by a Sersic profile (of the volume density with a total mass 
unity). Integrating the density, we find the central surface density ^(O)/^ 1 ~ 2goM/G ~ 
1011 ~ 13OOM0pc~ 2 for these systems, independent of the initial gas dispersion a and the 
feedback parameter F. Eventually /* = 1 when the gas is exhausted by star formation, and we 
form a stellar system with a bulge-like density profile. Our model resembles real bulges in term 
of the central brightness 7*(0) and the Faber- Jackson- like relation ~ at between the total 
mass and stellar dispersion. 

During the star burst (SB), the central Black Hole (BH) accretes, the SB and BH emit photons 
with a luminosity Lsb+bh, which diffuse out of the gas sphere while keeping the gas isothermal. 
The momentum deposite rate — SB + BH Q f photons drives an overall feedback force acting on 
the gas, which can be computed by 

2L SB +BH _ „,.._ f°° ... 2w 2s r 10(l-/.)r«7j 



F(t) 



y°° d(4^r 2 )(pa 2 )r ~ 10( | + p )r ^ ~ 2(1 - f.)M?g DM (4.5) 



insensitive to F if T ^ 1. The luminosity and force are maximum initially and die out as the 
star burst finishes. If the maximum luminosity at the onset of star burst is contributed fifty-fifty 
between the luminosity of the SB and the Eddington luminosity of the BH, then we obtain 



Mbh L.sl 



10 s M Q 10 13 L© V200kms 



which agree with the observed scaling relations of the stellar dispersion cr* with the BH mass 
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Figure 1. shows mass distribution of an isothermal gas sphere embedded in a NFW dark 
matter cusp of a uniform gravity gnu ~ 10~ 10 m/sec 2 for models (from bottom to top) with 
increasing dimensionless central gas pressure po = 1,3,10 (in solid green), and po = 30,100 
(black dashed). Panel (a) shows the models in logp(m) vs logx(m) (rescaled density vs. rescaled 

radius xim) — .'' ). Note how the black curves are above the green curves at small radii, 

but dip below the green curves at large radii, a feature of reaching a maximum in total gas mass 
at the critical pressure (po = 30). Also shown is a Sersic (n = 1.2) profile (red circles). Panel (b) 

shows lnp(m) (the logarithm of rescaled gas density or pressure, p(m) = ^ A f?f „, ) as function 

9 DM > 

of the rescaled enclosed gas mass m — M/(a 4 /G/qdm), cf. eq. 3.4). Note the total gas mass 
m increases with po until the critical value po ~ 30, afterwards the mass decreases with central 
pressure. 



and the star burst luminosity respectively. Assuming the usual SB efficiency of 0.001, the star 
formation time scale is O.OOIM^ c 2 /Lsb ~ O.OAHq 1 , comparable to the free-fall time scale 
~ O.lGyr. The short time scales and high luminosity are consistent with the assumption of 
violent feedback f>l. 
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